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AN EXPERIMENTAL INVESTIGATION 
OF THE MIXING OF COMPRESSIBLE -AIR JETS 
IN A COAXIAL CONFIGURATION 

By James M. Eggers and Marvin G. Torrence 
Langley Research Center 

SUMMARY 

An experimental investigation has been made to obtain data on the mixing of paral- 
lel, circular, compressible-air jets in a coaxial configuration. A Mach 1.3 outer jet and 
a near -sonic central jet, both exhausting to a quiescent atmosphere, were studied. Data 
were also obtained with both jets operating at subsonic Mach numbers. Radial distribu- 
tions of total pressure, static pressure, and mass-fraction concentration (air in air) of 
the central- jet gas were determined at downstream axial stations both in the near -jet and 
far -jet flow fields. Concentration measurements were made by using a tracer -gas tech- 
nique. The central-jet gas was a 1-percent mixture of ethylene in air. At each test con- 
dition, the total temperature of both jets was approximately 300° K (80° F). The Reynolds 
numbers for the supersonic and subsonic outer-jet flow were approximately 42.1 x 10^ and 
13.9 X 1()6 per meter (12.8 x 10® and 4.24 X 10^ per foot), respectively. 

Results of the investigation indicate that good correlation with the experimental air 
mixing data can be obtained by using an eddy viscosity model based on a representative 
width of the mixing zone and on the center -line mass flow per unit area. The concentra- 
tion decay in the region well downstream of the concentration potential core was found to 
be inversely proportional to the square of the axial coordinate. The accuracy of the con- 
centration measurements was found to be significantly affected by sampling technique and 
sampling-probe design. 


INTRODUCTION 

An important application of the turbulent mixing of compressible parallel jets is in 
the design of hypersonic ramjet-engine combustors. Although injection of the fuel per- 
pendicular to the airstream offers the advantage of greater penetration, the downstream 
momentum of parallel-injected fuel can contribute significantly to the engine thrust. (See 
ref. 1.) Since both the fuel distribution and combustor length will have an effect on the 
ramjet-engine performance, development of analytical techniques to predict the combustor 



flow field in order to optimize the combustor design is desirable. Current analytical 
techniques employ eddy viscosity formulations derived from nonreacting flow data, since 
extensive measurements representative of the severe environment of the ramjet combus- 
tor are not available. The assumption is made in the analysis that the mixing of the 
reacting flow field is adequately expressed by the parameters (for example, velocity, 
density, and mixing -zone width) in the nonreacting -flow eddy viscosity formulation. 

An analysis of the mixing of the turbulent flow fields has recently been performed 
by employing the laminar -flow boundary -layer equations modified by replacing the lam- 
inar viscosity with the eddy viscosity and by replacing the laminar Prandtl, Lewis, and 
Schmidt numbers by their turbulent counterparts. In reality, the eddy viscosity is a com- 
plex function of the local properties in the flow field, but an assumption of uniform eddy 
viscosity across the mixing zone is usually applied to facilitate computations. Since the 
introduction of Prandtl's mixing length theory in 1925 and his later simplified eddy vis- 
cosity theory in 1942, which was based on the difference in velocity across the mixing 
zone (see ref. 2), many eddy viscosity formulations have been proposed (see refs. 3, 4, 5, 
6, 7, and 8). No universally accepted eddy viscosity expression is currently available. 
Recent attempts have been made to gain sufficient knowledge of the eddy viscosity distri- 
bution to develop improved eddy viscosity models and also to determine the associated 
turbulent Prandtl, Lewis, and Schmidt numbers by calculating these parameters directly 
from experimental data. (See refs. 1, 9, 10, and 11.) Extreme difficulties, however, 
arise in determining these parameters by this direct method, since derivatives must be 
obtained from the experimental data. The result is a significant scatter of the values 
obtained; thus this method was not used herein. 

The current experimental investigation was designed to supply data to aid in devel- 
oping a turbulent mixing analysis and the associated eddy viscosity formulation. Com- 
pressible, turbulent mixing data were obtained from coaxial, parallel air jets which 
exhausted to a quiescent atmosphere. A Mach 1.3 outer jet, approximately 17.8 cm 
(7.0 in.) in diameter, surrounded a near-sonic central jet with a 2.443-cm (0.962-in.) 
diameter. In a later stage of the experiment, this central jet was replaced by a smaller 
1.036 -cm -diameter (0.408-in.) central jet to permit data acquisition at far downstream 
locations in terms of central-jet diameter in the flow field. Data were also obtained with 
the smaller central jet exhausting at Mach 0.47 and the outer jet at Mach 0.58. In all 
tests, the total temperature of the central and outer jet was approximately 300° K (80° F). 
The Reynolds numbers for the supersonic and subsonic outer -jet flow were approximately 
42.1 x 10® and 13.9 X 10® per meter (12.8 X 10® and 4.24 x 10® per foot), respectively. 
Radial distributions of total pressure, static pressure, and concentration of the central- 
jet gas (the latter by the tracer -gas technique) were obtained at various locations both in 
the near -jet and far -jet flow fields. 
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Data correlation was performed by employing the analysis of reference 7 in con- 
junction with an improved eddy viscosity model. This analysis was selected to correlate 
the data because of its immediate applicability to combustion problems. A discussion of 
the inaccuracies of gas samples taken from flowing streams is presented in the appendix. 

SYMBOLS 

The units used for the physical quantities defined in this paper are given both in the 
International System of Units (SI) and in the U.S. Customary Units. Factors relating the 
two systems are given in reference 12. 

A area 

b a mixing -zone width defined as the radial distance between the points where 

the velocities are uj and U 2 

D central-jet inside diameter (see fig. 1) 

I the ratio of integrated central-jet flow rate calculated by using the measured 

pressures, temperatures, and mass concentrations, the latter obtained by 
use of the high -flow -rate sampling system, to the metered central-jet flow 
rate 

I' the ratio of integrated central-jet flow rate calculated by using the measured 

pressures, temperatures, and mass concentrations, the latter obtained by 
use of the low-flow-rate sampling system, to the metered central-jet flow 
rate 

k empirical constant employed in eddy viscosity formulations 

M Mach number 

m local mass flow rate 

m c jet flow rate calculated by using equation (1) 

N]_,e turbulent Lewis number, the product of the turbulent diffusion coefficient 

and the constant -pressure specific heat divided by the turbulent thermal- 
conductivity coefficient 
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Npp turbulent Prandtl number, the product of the constant -pressure specific heat 

and the turbulent viscosity divided by the turbulent thermal -conductivity 
coefficient 

Ng c turbulent Schmidt number, the ratio of the turbulent Prandtl number to the 

turbulent Lewis number 

p static pressure 

total temperature 

u velocity 

uj velocity defined by equation (5) 

U 2 velocity defined by equation (6) 

u x free -stream velocity modified for static -pressure changes at the x location 

(see eq. (3)) 

V local velocity modified for static pressure and free -stream velocity deviations 

(see eq. (3)) 

Vp velocity modified for static -pressure changes (see eq. (3)) 

x axial coordinate (see fig. 4) 

Xj potential core length based on concentration 

y radial coordinate (see fig. 4) 

z a mixing -zone width based on the distance from the jet center line to the 

point where the velocity is the mean velocity in the mixing zone 

a local mass concentration obtained by using the high-flow-rate sampling 

system 

a;' local mass concentration obtained by using the low -flow-rate sampling 

system 


4 



eddy viscosity 


e t 

ef- 0 eddy viscosity evaluated at the jet exit station 

p density 

Subscripts: 

a evaluated in the quiescent atmosphere 

CL evaluated on the jet center line 

e evaluated in the outer -jet flow in the region of uniform velocity at the jet 

exit station 

j evaluated at the central -jet exit station and on the center line 

max maximum 

oo evaluated in the outer -jet flow external to the mixing zone in the region of 

uniform velocity 


APPARATUS 

Nozzle 

A sketch of the nozzle section is shown in figure 1. The nozzle contour was 
designed by using a modification of the technique of reference 13. Reference 13 derives 
analytical expressions for the design of axisymmetric Laval nozzles and expresses the 
desired contour as a function of the radial and axial coordinates. The analytical expres- 
sions of reference 13 were modified to predict the area ratio as a function of axial dis- 
tance for the desired Mach 1.3 plug nozzle. This design technique gave a reproducible 
method of generating a contour and resulted in a nozzle which produced an acceptable 
exit -velocity profile. However, significant shock waves are evident in the schlieren 
photographs of figure 2. Pressure measurements indicated the central-nozzle lip shock 
corresponded to about a 2.5° turning of the outer -jet flow, or approximately a 6 -percent 
change in velocity of the outer jet. Similarly, the outer -nozzle lip shock was found to 
correspond to about a 1.5° to 2° turning of the outer-jet flow. The outer -nozzle lip shock 
caused difficulties in data reduction as it impinged upon the central-nozzle flow and 
imposed a pressure rise upon the center region of the jet flow field. Further discussion 
of the effects of the shocks will be found in the section entitled "Presentation of Data." 
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Test Procedure 


The nozzle section of figure 1 was connected to a blower system with sufficient 
capacity to permit continuous operation of the nozzles. The nozzle exhaust pitot pres- 
sures were computed from the design exit Mach numbers by assuming that the exit static 
pressure was atmospheric. Reference pressures were monitored on mercury manom- 
eters to insure that the computed pitot pressures were maintained. Ethylene -tracer flow 
was metered by a variable-area orifice meter and injected into the central-jet air flow 
far upstream of the nozzle section to insure uniform concentration at the nozzle exit. 

The ethylene -tracer flow rate was approximately 1 percent by mass of the central- jet 
air -flow rate. 

After steady-state conditions were attained, point by point traverses of the flow field 
with a rake containing a pitot probe and static -pressure probe (one probe trailed the other 
along the line of the traverse) were made throughout the unconfined flow field. Pertinent 
information on the rake and probes is shown in figure 3. A flow-field schematic with the 
nomenclature used is shown in figure 4. A single -pass schlieren system was used to 
visually observe and photographically record the jet flow fields. Other supporting mea- 
surements recorded include probe location, total temperature of both jets, static pressure 
from the taps around the inside of the outer-nozzle lip (see fig. 1), barometric readings, 
and central-jet flow rate. The central-jet air flow was measured by a calibrated thin- 
plate -orifice meter. 

Gas samples were extracted from the flow field through the pitot probe by using only 
the pitot pressure as the pumping force. Approximately 3 percent of the sample flow was 
routed through a gas chromatograph and analyzed for ethylene content. The remainder of 
the sample flow was bypassed and discharged to the atmosphere. (Further details on 
sampling technique are discussed in the appendix.) Both the flow through the gas chroma- 
tograph and the bypass flow were measured by flow meters. Data-acquisition time was 
approximately 3 minutes per data point and was essentially fixed by the time of the gas- 
chromatograph-analysis cycle. A discussion of concentration measurement technique 
and theory by use of the gas chromatograph may be found in reference 14. Data were 
obtained for the three test conditions shown in the following table: 


Test 

condition 

Central -jet 
diameter 

Mj 

M e 

U j 

u e 

T tJ 

T t,e 


1 

2.443 cm (0.962 in.) 

0.942 

1.302 

0.736 

1.0 

49 

2 

1.036 cm ( .408 in.) 

.813 

1.268 

.648 

1.0 

114 

3 

1.036 cm ( .408 in.) 

.47 

.58 

.827 

1.0 

85 
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Several samples of the central-jet gas were extracted from the piping upstream of 
the nozzle exit at the beginning and end of each test. The chromatograph output from the 
analysis of these samples was set to correspond to full-scale pen deflection on a recorder. 
This procedure established a central-jet reference concentration based on ethylene which 
corresponded to 100 percent mass concentration of central -jet gas. (Although the gas 
chromatograph analysis yields mole concentration, for gases of nearly equal molecular 
weights as between ethylene and air, mole and mass concentration are equivalent.) The 
chromatograph output was experimentally verified to be linear with concentration. There- 
fore, a sample which produced a 50 percent full-scale pen deflection from the chromato- 
graph output contained 50 percent as much ethylene as the central-jet gas, and thus the 
mass concentration of the sample was 50 percent of that of the central-jet gas. The 
reduced -concentration data are reported in terms of local mass fraction a of central- 
jet gas. For clarity, the units of a may be thought of as the ratio of pounds central- 
jet gas per pound of mixture. 


Accuracy of Results 

Based on the chart resolution of the recorded pitot pressures and an uncertainty of 
the total-temperature measurement of ±6° K (±10.8° F), it is estimated that the velocity 
data are accurate to ±2 percent. This accuracy does not include any uncertainty due to 
the static pressure measurements. In the data of test conditions 1 and 2, the static pres- 
sure was not uniform in either the axial or radial direction. 

One means of assessing the accuracy of the data is to apply the principle of conser- 
vation of mass to the central-jet flow. The following equation is evaluated for each axial 
station at which data are obtained for this assessment: 


m, 


= ^ Q!(m/A)dA 


( 1 ) 


where m/A is the local mass flow per unit area evaluated from the total temperature, 
pitot pressure, and static pressures measured, a is the local mass concentration, and 
A is the area of the flow field over which a is nonzero. If the data are correct, m c 
will be equal to the central-jet flow rate. Differences between m c and the experimen- 
tally determined central-jet flow rate obviously are due to experimental error. A com- 
parison between these two mass -flow quantities is an essential criterion in assessing the 
accuracy of concentration measurements because of uncertainties in obtaining represen- 
tative samples from flowing streams. 

The application of equation (1) to the data presented herein is reported in the appen- 
dix. Large errors in concentration measurements taken from binary streams may occur, 
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the probe design and sampling technique employed having a significant effect upon the 
results. The actual physical mechanism which causes unrepresentative sample collection 
is not known, but satisfactory measurements were obtained with the probe and sampling 
technique used in this investigation. The results suggest that the erroneous concentration 
measurements are related to the local turbulence level in the flow field. 

Theory 

The analysis of reference 7 employs equilibrium chemistry, transformation, and 
explicit finite difference techniques to compute turbulent mixing and combustion of paral- 
lel unconfined streams of hydrogen and air. Both two-dimensional and axisymmetric con- 
figurations are treated, the method differing only slightly in the transformation employed. 
The analysis has provision for axial pressure gradients and nonunity turbulent Prandtl 
and Lewis numbers. Principal inputs to the mixing program are distributions of velocity, 
static temperature, and mass concentration corresponding to equal increments of the 
stream function. This input requirement is inconvenient for other than step profiles as 
considerable calculation is required to determine the stream -function distribution. The 
analysis initially contained an eddy viscosity formulation attributed by reference 15 to 
Zakkay. The Zakkay formulation was based on experimental concentration data obtained 
downstream of the potential core and was intended to be used only in far -field calculations. 
It was emphasized that the Zakkay eddy viscosity formulation predicts a concentration 
decay proportional to l/x^. Zakkay' s eddy viscosity formulation has the form: 

e t = kz(pu) CL (2) 

In equation (2), ej- is the eddy viscosity, z is the mixing -zone width based on the 
distance from the jet center line to the point where the velocity is the mean velocity in 
the mixing zone, p is the density, and u is the velocity. A value of 0.04 for the empir- 
ical constant k was suggested in reference 7 to be used in conjunction with turbulent 
Prandtl and Lewis numbers of unity, step velocity -profile input, and mixing of hydrogen 
and air. More detailed information on the mixing program is contained in references 7 
and 15. 

In order to analyze the air -in-air mixing data obtained herein, the central jet was 
input as standard air and the outer jet was input as 100 percent nitrogen. The oxygen 
concentration could be interpreted as a tracer gas in the same manner that ethylene was 
used in obtaining the experimental data. 

PRESENTATION OF DATA 

The experimental data of test conditions 1 and 2 are presented in figures 5 and 6, 
respectively. It was observed during reduction of the data for test conditions 1 and 2 that 
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the experimental static -pressure distributions were not uniform in either the axial or 
radial direction. A nonuniform, radial static -pressure distribution is not provided for 
in current mixing analyses; therefore, it was necessary to devise a data-reduction method 
whereby the nonuniformity of static pressure could be eliminated from the data. 

Two possible methods to reduce the data were suggested. The first method was to 
assume uniform static pressure equal to atmospheric pressure throughout the flow field. 
The second method was to use the static-pressure distribution from the nozzle exit to 
reduce the data at the downstream stations. The second method was chosen because 
retention of a true representation of the initial mixing conditions was desirable and an 
assumption of uniform static pressure was unrealistic. Center -line static -pressure 
ratios, that is, measured static pressure divided by atmospheric pressure, for test con- 
dition 1 are tabulated as follows: 


x/D 

0 

11.0 

17.2 

25.0 

30.0 

36.0 

49.0 

P/Pa 

1.11 

1.18 

1.15 

1.15 

1.15 

1.11 

1.04 


The static pressure at all axial locations except the jet exit generally decreased with 
increasing radial dimension until ambient-pressure levels were attained. 

The outer -jet -exit static pressure was determined from the static taps around the 
outer -nozzle exit and from the probe static -pressure measurement. Where shock waves 
did not intercept the static -pressure probe, the probe and static -tap measurements were 
in good agreement. The outer-jet-exit static -pressure ratio P/P a was 0.98 for test 
condition 1. 

The center -line velocity distribution for test condition 1 obtained by employing the 
measured static pressure is labeled experiment in figure 5(a). Similar distributions of 
static pressure were observed in the data of test condition 2. However, even more dis- 
tinct effects of compressions and expansions were evident and are reflected in the more 
sizable scatter of the experimental center-line velocity data of figure 6(a). These data 
were computed by using the measured static pressure. As previously stated, it was 
desirable to retain the measured nozzle exhaust conditions and reduce all downstream 
profiles to the static -pressure level at the nozzle exit. As no precise method existed for 
incorporating this change, because of the radial as well as the axial static -pressure vari- 
ation, the somewhat arbitrary decision was made to make the static -pressure change along 
streamlines. Further analyses indicated that for the given conditions the streamlines, 
actually stream -tube surfaces, were reasonably parallel over the data range. Thus, the 
data-reduction method allowed static -pressure changes to be made along lines of constant 
radius. Therefore, in all data except those at the jet exit station, the local measured pitot 
pressure was used in conjunction with the static pressure measured at the jet exit at the 
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same radial location to compute local Mach numbers. Velocities were then computed by- 
using an average total temperature of 300° K (80° F) and a specific -heat ratio of 1.4. 

In addition to the pressure modification, a second modification was incorporated into 
the velocity data which compensated for variation of the free -stream velocity. These 
variations were due to experimental error, to viscous and shock losses and, at the most 
downstream stations, to penetration of the outer mixing boundaries of the outer jet and 
surrounding air into the mixing zone (see fig. 4). Similarity of profiles was applied 
through the following equation: 


V - uj _ V p - uj 
u e - Uj u x - Uj 


( 3 ) 


where V is the new local velocity modified for pressure and free-stream velocity devi- 
ation, Uj is the center -line jet exit velocity, u e is the velocity evaluated in the outer - 
jet flow at the jet exit in the region of uniform flow, u x is the free-stream velocity mod- 
ified for static -pressure changes at the x location, and Vp is the local velocity modified 
for static -pressure changes. All velocity data for test conditions 1 and 2, except the basic 
data points in figures 5(a) and 6(a), have been modified for both static-pressure changes 
and free-stream velocity deviations according to the previously discussed methods. 

The resulting velocity profiles corresponding to test condition 1 are shown in fig- 
ures 5(c) and 5(d). The magnitude of the free-stream velocity modification is less than 
2 percent at x/D less than 30, less than 3 percent at x/D of 30 and 36, and less than 
10 percent at x/D of 49. The uncertainty of the overall pressure and velocity modifi- 
cation is reflected in scatter in the velocity data downstream of the inner nozzle lip (for 
y/D near 0.5) as is evident in the data of figures 5(c) and 5(d). 

The resulting velocity profiles corresponding to test condition 2 are shown in fig- 
ures 6(c), 6(d), 6(e), and 6(f). Pressure and free-stream velocity modifications were 
made as in the data of test condition 1. For test condition 2 data, the magnitude of the 
free-stream velocity correction was less than 2 percent at x/D less than 102.8 and less 
than 4 percent at x/D of 102.8 and 114.2. As in the test condition 1 data, uncertainties 
of the modifications are evident in the velocity profiles for values of y/D near 0.5. 

Concentration distributions for test condition 1 are presented in figures 5(b), 5(e), 
and 5(f). Concentration distributions for test condition 2 are presented in figures 6(b), 

6(g), 6(h), 6(i), and 6(j). The concentration data for both test conditions were consistent 
except for scatter in the data of test condition 1 at an x/D of 21.4. The four concentra- 
tion values at an x/D of 21.4 in figure 5(b) represent data obtained from different runs 
and on different days. No satisfactory explanation for the scatter was found. 
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Experimental data corresponding to test condition 3 are shown in figure 7. The 
center -line velocity distribution and the velocity profiles for test condition 3 are shown 
in figures 7(a) and 7(c) and (d), respectively. No modification for pressure or free- 
stream velocity variations was applied to the test condition 3 data. Local measurements 
of pitot and static pressures were used to obtain Mach number distributions. A mean 
total temper ature was then employed to obtain velocity distributions by using a specific - 
heat ratio of 1.4. 

The experimental center -line concentration distribution and the experimental 
concentration profiles for test condition 3 are shown in figures 7(b) and 7(e) and (f), 
respectively. 


Data Analysis 

As previously discussed theoretical analyses of the experimental flow fields were 
performed by employing the mixing analysis of reference 7 and the associated computer 
program described in reference 15. As may be seen in figures 5(c), 6(c), and 7(c), the 
jet exit flow consisted of a large percentage of boundary -layer flow. Therefore, the near- 
field velocity profiles are not similar, the initial mixing is wakelike in nature, and a step- 
function velocity -profile input to the mixing program could not be expected to satisfactorily 
predict the mixing in the near field. Thus, it was necessary to employ a representation of 
the experimental velocity profiles as inputs to the mixing program. The solid line in fig- 
ures 5(c), 6(c), and 7(c) at an x/D of zero is representative of the input velocity profiles. 
Between 23 and 39 input points were employed to represent the profiles, the exact number 
depending upon the particular test -condition data being analyzed. Since a velocity of zero 
at the inner -nozzle lip (y/D = 0.5) violates the mathematics of the analysis, a low but finite 
velocity value was used as an input at this location. Step profiles of concentration were 
employed as inputs for all computations performed, as shown in figures 5(e), 6(g), and 7(e). 

Initial computations were performed by using the constants recommended in refer- 
ence 7, that is, k = 0.04, Np r = 1.0, Nl 6 = 1.0, and the Zakkay viscosity model (see 
eq. (2)). The results of these calculations are shown as dashed lines in figures 5(a) and 
5(b), 6(a) and 6(b), and 7(a) and 7(b) for test conditions 1, 2, and 3, respectively. It is evi- 
dent that the Zakkay viscosity model in conjunction with k = 0.04, greatly overpredicts 
the rate of mixing of the air jets for all three test conditions. The theoretical solution 
for other values of k may be generated from the solution for a given k value by scaling 
x/D by the ratio of the k's involved. A more nearly correct value of k to be used in 
conjunction with the Zakkay viscosity model would be on the order of 0.004 for the far- 
field velocity data. Initial calculations indicated that satisfactory correlation of the near- 
field mixing data of test condition 1 could be obtained by reducing the constant employed 
with the Zakkay viscosity model to 0.0025. However, as the Zakkay viscosity model was 
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formulated for far -field calculations and contains a mixing-zone width which is represen- 
tative of the mixing zone only in the far field, it is not surprising that continuous calcula- 
tions throughout the flow field were unsatisfactory. The center -line velocity distribution 
predicted by the Zakkay viscosity model and by using k = 0.0025 is shown in figure 5(a) 
for the data of test condition 1. It is evident that the Zakkay model does not produce a 
rapid enough increase in eddy viscosity at the end of the near field to correctly predict 
the far -field data in one continuous calculation. 

Illustrative of the requirement for increasing eddy viscosity at the end of the near 
field are the preliminary correlations of the test condition 1 data which were presented 
in reference 16. (The preliminary data of ref. 16 cannot be compared in magnitude with 
those data presented herein, as a different data -reduction and correlation technique was 
used.) In reference 16, the Zakkay viscosity model was employed with a k of 0.0025 
in the near field (x/D £ 16.2) with a step change of k to 0.006 at the x/D = 16.2 loca- 
tion. Thus, the eddy viscosity was increased 2.4 times in the region near the end of the 
velocity potential core. Good correlation with the experimental data was attained; how- 
ever, the step change in k with associated discontinuous gradients of concentration and 
velocity together with the difficulty in determining the proper location in the flow field to 
make the step change in k made the calculation technique unsatisfactory. 

The difficulties in correlating the data by using the step change in k made it evi- 
dent that a method which would permit continuous calculations throughout the flow field was 
highly desirable. As a result, the following eddy viscosity expression was formulated: 

e t = kb(pu) CL (4) 

In equation (4), is the eddy viscosity, k is the empirical constant, (pu)^ is the 
mass flow per unit area on the jet center line, and b is the mixing -zone width between 
the points where the local velocity is uj and U2 as given by equations (5) and (6). 


Ui = U^ + 0.95 (u cl - u^) 

(5) 

u 2 = u M + 0.50 (u cl - Uoo) 

(6) 


The mixing-zone width b thus is based on a width of the mixing zone rather than a dis- 
tance from the jet center line as in the Zakkay model, and, therefore, has physical signif- 
icance in both the near -jet and far -jet flow fields. 

Calculations for all three test conditions were performed by using the eddy viscosity 
expression given by . equation (4). The resulting theoretical center -line velocity distribu- 
tions are shown by the solid lines in figures 5(a), 6(a), and 7(a) for test conditions 1, 2, 
and 3, respectively. Good correlation was obtained for the data of all three test condi- 
tions; however, different values of the empirical constant were required for the data of 
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each test condition. The empirical constants required were 0.0098, 0.0078, and 0.011 
for the data of test conditions 1, 2, and 3, respectively. The theoretical velocity profiles 
are compared with the experimental velocity profiles for test condition 1 in figures 5(c) 
and 5(d), for test condition 2 in figures 6(c), 6(d), 6(e), and 6(f), and for test condition 3 
in figures 7(c) and 7(d). Good agreement between the theoretical velocity profiles and 
the experimental data is evident in the figures for the data of all three test conditions. 
Emphasis is placed on the center -line distributions, however, since a reasonable fit 
throughout the flow field indicates that only one well-selected experimental data point on 
the center line is necessary to determine the experimental constant k for given initial 
mixing conditions. It should also be pointed out that the theoretical velocity solutions 
for all three test conditions are independent of the magnitude of the turbulent Prandtl and 
turbulent Lewis numbers, as no significant transfer of heat is involved. 

The mass concentration of central -jet gas on the center line as predicted by using 
the analysis of reference 7 and the eddy viscosity model of equation (4) is presented in 
figures 5(b), 6(b), and 7(b) for test conditions 1, 2, and 3, respectively. Note that solu- 
tions for various turbulent Schmidt numbers Ng c (the turbulent Schmidt number is the 
ratio of the turbulent Prandtl number Np r to the turbulent Lewis number are pre- 

sented, as examples of the effect of the variation of this parameter are quite limited in 
the literature. It is evident that increasing the Schmidt number by a factor of four (from 
0.25 to 1.00) is equivalent to a reduction of the experimental constant by a factor of two 
to three depending upon its location in the flow field. However, a change in k also 
would cause a shift in the predicted velocity distribution. Thus, the correlation technique 
requires, first, employing the experimental center -line velocity distribution to determine 
the experimental constant k and then employing the center -line concentration distribu- 
tion to obtain the most appropriate value of turbulent Schmidt number. Whether the tur- 
bulent Prandtl number or turbulent Lewis number is 0.60 or 1.00 cannot be determined 
from the data herein. However, the turbulent Schmidt number is sufficient to determine 
the uniqueness of the concentration distribution along with the empirical constant k. 
Inspection of figures 5(b) and 6(b) indicates that reasonable correlation of the center -line 
concentration distribution is obtained by using a turbulent Schmidt number on the order 
of 0.60. However, for the test condition 3 data (fig. 7(b)), the turbulent Schmidt number 
required to correlate the data is somewhat higher. Theoretical calculations for the test 
condition 3 data were made by using a turbulent Schmidt number of unity. No reason 
for the differences in the turbulent Schmidt number required for correlation is evident. 
Inadequacies in the eddy viscosity expression, for example, too fast or too slow mixing 
at a location in the flow field, will affect not only the velocity but also the predicted con- 
centration distribution. Particularly evident near the end of the concentration potential 
core is the too rapid mixing of the eddy viscosity model (or possibly the uncertainty in 
turbulent Schmidt number) which causes underprediction of the center -line concentration. 
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Theoretical concentration profiles are compared with the experimental data in figures 5(e) 
and 5(f) for the data of test condition 1, in figures 6(g), 6(h), 6(i), and 6(j) for the data of 
test condition 2, and in figures 7(e) and 7(f) for the data of test condition 3. Good agree- 
ment between theory and the experimental data is evident in the figures for all three sets 
of experimental data. It is concluded that an eddy viscosity model based on a representa- 
tive width of the mixing zone and the center -line mass flux per unit area results in good 
correlation with the air -mixing data reported herein. However, the experimental con- 
stant required to correlate the data apparently depends upon the experimental conditions 
of the data considered. 

It has been reported in reference 11 that the decay of center -line mass concentra- 
tion of axisymmetric jets is inversely proportional to x^. This decay rate was based on 
experimental data of jets of various gases, including argon, carbon dioxide, helium, and 
hydrogen. Figure 8 indicates that the observed concentration decay in the region well 
downstream of the concentration potential core is approximated by the inverse of the x^ 
decay rate. In the data of test conditions 2 and 3, the use of the inverse of the x^ decay 
rate to estimate concentration values just downstream of the end of the concentration 
potential core will result in substantial error for mass concentrations of 50 percent and 
greater. It was also reported in reference 11 that the concentration potential core length 
may be computed to a first approximation by the following equation: 



(7) 


(Note the concentration potential core length x^ has been nondimensionalized by the jet 
diameter here rather than the jet radius as in ref. 11. This changes the value of the con- 
stant from 22 to 11 as shown in eq. (7).) This equation was derived by extrapolation of 
the concentration decay curve (a ~ l/x2) to 99 percent mass concentration. A comparison 
of the concentration potential core lengths predicted by equation (7), by the l/x2 extrap- 
olation, and by a smooth extrapolation, the latter of which appears to give the best repre- 
sentation of the data, is shown in the following table: 


Test 

Concentration potential core length, x/D 

condition 

Equation (7) 

l/x2 extrapolation 

Smooth extrapolation 

1 

9.4 

17.5 

13.5 

2 

8.9 

17.0 

10.0 

3 

10.0 

14.5 

6.0 


It is concluded that the concentration potential core length is not accurately predicted by 
assuming a variation proportional to the square root of the ratio of the central-jet mass 
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flow per unit area to the outer -jet mass flow per unit area. It is observed that the data 
from which equation (7) was formulated have a scatter of 30 percent and -50 percent rela- 
tive to values determined by equation (7), thus accurate predictions could not be expected. 
Concentration potential core lengths obtained by a smooth extrapolation are shown as 
dashed lines in figure 8. The experimental center -line jet exit Mach numbers were 0.942, 
0.813, and 0.47 for test conditions 1, 2, and 3, respectively. A trend of decreasing con- 
centration potential core length with decreasing central -jet Mach number is noted. This 
trend is consistent with the smooth extrapolations of figure 8 and with the known fact that 
lower Mach number streams mix faster and, thus, are expected to have shorter potential 
core lengths. 


Eddy Viscosity Distributions 

Since the success of any turbulent mixing analysis depends strongly upon the ade- 
quacy of the particular eddy viscosity expression employed, it is instructive to examine 
the eddy viscosity distribution obtained in the correlation of the turbulent mixing data pre- 
sented herein. In figure 9, the variation of the ratio of the eddy viscosity at each x/D 
to the eddy viscosity at the jet exit with x/D is presented. Similar trends of the eddy 
viscosity increasing at an increasing rate near the jets exit (x/D < 20) and of the eddy 
viscosity varying linearly with x in the far field are evident for the data of all three test 
conditions. It is apparent in figure 9 that in the far field for a given x/D the curve for 
test condition 1 is approximately 2.36, the nozzle -diameter ratio, times the curve for test 
condition 2. (Note that the central-jet mass flow per unit area was approximately equal 
for the data of test conditions 1 and 2.) Also, for a given x/D the curve for test condi- 
tion 2 is approximately 1.80, the ratio of mass flows per unit area, times the curve for 
test condition 3. (Note that the jet diameter was the same for test conditions 2 and 3.) 
Thus, in the far field, the eddy viscosity scales as pjUjD; however, because the initial 
velocity profiles are nonsimilar, the near field will not scale in a similar manner. Initial 
values of the eddy viscosity for the three test conditions are shown in figure 9. 

Although the eddy viscosity expression employed correlates the data presented sat- 
isfactorily, certain limitations are inherent in the formulation. For example, the mixing - 
zone width b is based on velocity difference; thus, as the velocity difference between 
two streams becomes very small, the mixing-zone width becomes indefinite. Further- 
more, whether the present eddy viscosity model may be satisfactorily applied to flow 
fields where large density gradients occur, such as in hydrogen-air mixing, is unknown. 
Although the empirical constant k required adjusting to correlate the data,- it is antic- 
ipated that application of this eddy viscosity model to a sufficient quantity of data will pro- 
duce a relationship which will enable calculation of k for given initial mixing conditions. 
Since the eddy viscosity is semiempirical in nature, it is too much to expect a single cor- 
relation to be valid for all conditions encountered. The eddy viscosity model employed 
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herein is similar to that employed in reference 17 in which an analysis of the mixing of 
unconfined jets with quiescent air was performed. Ip. addition the eddy viscosity model 
used herein was successfully employed in reference 18 to compute two-dimensional mixing 
of supersonic streams in the presence of shock waves. 

CONCLUSIONS 

As a result of the experimental and analytical study of the turbulent mixing of cir- 
cular, coaxial jets the following conclusions were obtained: 

1. An eddy viscosity model based on a representative width of the mixing zone and 
the center -line mass flow per unit area together with the experimentally determined con- 
stants resulted in good correlation with the air mixing data reported herein. 

2. Large errors in concentration measurements taken from binary streams may 
occur, the probe design and sampling technique employed having a significant effect upon 
the results. The actual physical mechanism which causes unrepresentative sample col- 
lection is not known, but satisfactory measurements were obtained with the probe and 
sampling technique used in this investigation. The results suggest that the erroneous 
concentration measurements are related to the local turbulence level in the flow field. 

3. The concentration decay in the region well downstream of the concentration poten- 
tial core is well approximated by an inverse of the square of the axial coordinate. 

4. The concentration potential core length is not accurately predicted by assuming 
a variation proportional to the square root of the ratio of central-jet to outer-jet mass 
flow per unit area. 

5. Turbulent Schmidt numbers required to correlate the concentration data were 
found to be approximately 0.6 for the high-velocity data and somewhat larger for the low- 
velocity data. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 1, 1969, 

722-03-00-04-23. 
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APPENDIX 


SAMPLING TECHNIQUE 

Large differences in concentration measurements may occur, the sampling tech- 
nique employed having a significant effect upon the values obtained. In this investigation 
the central-jet flow computed from the data by using equation (1) was found to be approx- 
imately 25-percent less than the metered central-jet flow. In reference 14 similarly 
evaluated accuracies of ±20 percent were found in mass flows when ethane was used as a 
tracer gas. In other investigations (refs. 19 and 20), an excess of the larger molecular 
weight gases were found in the gas samples. Recent literature (ref. 5) stresses the fact 
that the probe design can radically affect the concentration measurement. 

It was postulated that if the probe captured all flow in the stream tube in front of its 
sampling port, a minimum of sampling bias would occur. Up to this time only sufficient 
flow to supply the chromatograph (40 cm^/min) was allowed to pass through the probe and 
the remaining 97 percent of full stream -tube capture flow was spilled around the probe. 

The sampling system was modified to remove all possible restrictions from the sampling 
lines and an internally expanded sampling probe, shown in figure 3, was installed. The 
modified sampling system passed 40 cm3 /min through the chromatograph and bypassed 
approximately 1300 cm^/min during a checkout test. This total probe flow (1340 cm3/min) 
accounted for all but about 100 cm^/min of the computed full-capture probe flow for the 
checkout condition. New concentration readings were found to be higher than the pre- 
viously obtained data by approximately the 25 percent required to satisfy continuity of 
mass of the central-jet flow. 

Two sets of concentration data were acquired at test condition 2 by using the pitot 
probe shown in figure 3. The first set of concentration data was taken with approximately 
40 cm^/min flow through the sampling probe, and the corresponding local mass concen- 
tration of central-jet gas will be called cn'. The second set of data employed the modi- 
fied high -flow -rate sampling system discussed previously and the resulting local mass 
concentration will be called a. Switching from one sampling flow rate to the other was 
accomplished by opening a bypass valve, thus all other details of the two sampling tech- 
niques were identical. The ratio of the integrated to the metered central-jet flow, as 
determined by using the concentration a 1 , is given the symbol I'. Similarly, the ratio 
of the integrated to the measured central-jet flow, as determined by using the concentra- 
tion a, is given the symbol I. Values for these parameters at different values for x/D 
for the three test conditions are presented in table I. Values for I* are presented only 
for test condition 2. 
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TABLE I.- MASS- FLOW INTEGRATIONS 


Test condition 1 

Test condition 2 

Test condition 3 

x/D 

I 

x/D 

I 

I’ 

x/D 

I 

0 

0.96 

0 

0.94 


0 

1.02 

11.0 

1.05 

7.34 

1.01 

0.79 

7.34 

1.00 

17.2 

1.04 

13.5 

1.02 

.70 

13.5 

.96 

25.0 

.99 

20.8 

1.03 

.53 

20.8 

.99 

30.0 

1.05 

30.6 

1.00 

.52 

26.2 

.95 

36.0 

1.02 

36.7 

.99 

.45 

42.5 

1.00 

49.0 

1.03 

44.1 

1.05 

.44 

55.5 

1.03 



55.4 

.95 

.34 

73.5 

.94 



61.2 

.89 

.32 

85.0 

.83 



73.4 

.94 

.25 





84.9 

.87 

.18 





97.9 

.80 

.27 





102.8 

.79 

— 





114.2 

.88 

— 




It is evident from inspecting the values in table I that the low-flow-rate sampling 
technique accounted for as little as 18 percent of the total central-jet flow at x/D = 84.9. 
The variation of a.' /a with local Mach number is shown in figure 10(a). A ratio of 
a' /a of unity corresponds to no sampling-rate effect, and a value of zero indicates that 
the high -flow -rate sampling technique detected central-jet gas whereas the low -flow-rate 
sampling technique found none. The ratio a '/a may also be interpreted as an indica- 
tion of sampling error by using the low-flow-rate sampling technique. In figure 10(b) the 
ratio of a ' /a is plotted against the radial coordinate to show the location of each point 
in the flow field. The following observations are made with regard to the data of fig- 
ures 10(a) and 10(b): 

(1) The sampling error is not a function of Mach number alone, as different values 
of the sampling error a' /a are obtained at identical Mach numbers. 

(2) Errors occur in the subsonic flow as well as the supersonic flow. Thus, the 
explanation of reference 21, in which similar erroneous sampling results are attributed 
to a detached normal shock, is not applicable here. Similar sampling problems in sub- 
sonic flow were found in reference 5. 

(3) The sampling error is insensitive to local gradients in the flow field. At the 
x/D = 97.90 station, where the Mach number gradient and the concentration gradient are 
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small, the ratio a 1 /a varies from 0.125 to 0.625. Scatter at the far right of figure 10(b) 
is due to a' /a consisting of the ratio of two experimental quantities which are subject 
to sizable experimental error as their magnitudes become small. 

(4) The few points at x/D of 7.34 and 13.50, which indicate the least sampling-rate 
effect, are located near the center line of the jets and in the region of the potential core. 
(See fig. 10(b).) Thus, these points are in a region of relatively low turbulence. The 
fact that they have very small errors suggests that the sampling errors are related to the 
local turbulence level in the flow field. A similar effect of no sampling error near the 
jet exit was found in references 20 and 21. Although the data of reference 21 were gen- 
erated in a realm of extremely low pressure, the sampling effects are so strikingly sim- 
ilar to those found herein that the same mechanism is suggested. 

(5) A separation or pressure diffusion model, similar to those discussed in refer- 
ences 21, 22, and 23, in which light gas molecules are deflected due to streamline curva- 
ture near the sampling probe, and heavy gas molecules continue on into the probe, does 
not appear consistent with the observed fact that no erroneous samples are obtained within 
the potential core region. In particular, the data of reference 20 also show that in the 
region of the jet center line, the erroneous sampling results start at approximately the 
point where the potential core ends. 

It is concluded that large errors in concentration measurements taken from binary 
gas streams may occur and that the sampling technique employed has a significant effect 
upon the values obtained. It is also concluded that the actual physical mechanism which 
causes the sampling probe to obtain unrepresentative samples of gaseous mixtures is not 
known; however, the results suggest that the erroneous concentration measurements are 
related to the local turbulence level in the flow field. 
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Figure 1.- Sketch of nozzle test section. 





Central- nozzle 
shock (typical) 


(a) Knife edge normal to flow; D = 1.036 cm (0.408 in.); 
Mj = 0.81; M e = 1.27. 



(b) Knife edge normal to flow; D = 2.443 cm (0.962 in.); 
Mj = 0.94; M e = 1.30. 



(c) Knife edge parallel to flow-, D = 2.443 cm (0.962 in.); Id) Knife edge parallel to flow, time exposure; D - 2.443 cm (0.962 in.); 
Mj = 0.94; M e = 1.30. Mj = 0- 94 '. M e = L3a 


Figure 2.- Schlieren photographs of flow fields. 
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Figure 4.- Flow-field schematic. 
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(a) Test condition 1. 
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